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Abstract 

We compute the mutual information between finite intervals in two non-compact 
2d CFTs in the thermofield double formulation after one of them has been locally 
perturbed by a primary operator at some time ti^ in the large c limit. We determine 
the time scale, called the scrambling time, at which the mutual information vanishes 
and the original entanglement between the thermofield double gets destroyed by the 
perturbation. We provide a holographic description in terms of a free falling particle 
in the eternal BTZ black hole that exactly matches our CFT calculations. Our results 
hold for any time t^j. In particular, when the latter is large, they reproduce the bulk 
shock-wave propagation along the BTZ horizon description. 
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1 Introduction and Summary 


Consider a (l+l)d CFT in some state having non-trivial qnantnm correlations. If the system 
is pertnrbed at some instant of time and evolve nnitarily afterwards, it is natnral to ask 
whether there exists any time scale when its snbsystems become nncorrelated. Since the 
mntnal information Ia-b = Sa+Sb — Saub between two snch snbsystems A and B provides an 
npper bonnd for the connected two-point fnnctions of operators acting on these snbsystem^ 

[I] 

, {{OaOb) - {OA){OB)f . . 

- 2||C>aPI|Ob|P ’ 

it is natnral to stndy the vanishing of this qnantity to answer this qnestion. The stndy of 
the time dependence in this measnre of entanglement can help ns to nnderstand the time 
scales controlling how qnantnm systems get thermalized, which is one of the most important 
problems in non-eqnilibrinm physics. 

A particnlar sitnation of the above scenario is when a pertnrbation acts on a thermal 
state. In holographic theories, thermal states are believed to have a gravity dnal in terms 
of black holes PE]. Black hole physics snggests that the speed at which the system forgets 
initial conditions, i.e. the pertnrbation, is the fastest among all physical diffnsive processes. 
This gave rise to the notion of fast scramblers and the scrambling conjectnre BE]- The main 
goal of this paper is to provide a hrst principle derivation for the time scale at which this 
phenomenon occnrs for 2d CFTs in the large c limit in a concrete setnp which allows both 
CFT and holographic compntations. 

Recently, Shenker and Stanford considered an excellent and tractable setnp to stndy the 
fast scrambling phenomena in the context of an eternal black hole BEli- This involves 
a pair of non-interacting CFTs in an entangled state, the thermoheld donble state. Tracing 
any entire CFT Hilbert space, gives rise to a thermal density matrix in the remaining CFT. 
The pertnrbation is described by some bonndary CFT operator and its gravity dnal involved 
a shock-wave propagating in the black hole backgronnd. No matter how small the bonndary 
pertnrbation is, the bine shift of energies when this pertnrbation reaches the horizon snggests 
the existence of a non-trivial backreaction. 

In this work, we stndy snch a setnp for a pertnrbation localized in a point-like region, 
triggered by a primary operator in a given CFT. To obtain analytical resnlts, we consider 
2d large c CFTs and their gravity dnals given by a pertnrbation of the BTZ black hole [5]. 
Recent developments in the calcnlation of 4-pt fnnctions involving heavy and light operators 
in the large c limit of the dnal 2d CFT mm (see also PEIOE]) allow us to analytically test 
these ideas. 

Computations of time evolutions of entanglement entropy after local perturbation^ by 

^This bound is proved for finite dimensional Hilbert spaces. We are not aware of an extension of this 
result to QFTs/CFTs, but we expect it to hold when regulating and normalising appropriately the relevant 
quantities in the continuum limit. 

^Note that this setup looks similar to the local quenches in CFTs m- However, in the latter the local 


2 





primary operators have been formulated in [ElE] and have been applied to many examples 
for CFTs at zero temperature in [121 [HI IHl Ell E21 E3]- Entanglement entropy and mutual 
information at hnite temperature CFTs has been analyzed for integral CFTs in [21]. On the 
other hand, the holographic calculations of time-evolutions of entanglement entropy after 
local perturbations have been analyzed in [TKl EH ES] at zero temperature and in [21] at 
hnite temperature. In this paper we will extend the discussion of local excitations to the 
thermoheld double formalism of hnite temperature CFTs in the large c limit. 


Summary of results : Specihcally, we perturb the thermoheld double (TFD) state by a 
local primary operator -0 at time in the past and compute the mutual information between 
regions A and B belonging to opposite boundaries. We denote the two boundary times in 
the thermoheld double by ti = and t^ = t^. When measuring the mutual information at 
= t_i_ = 0, we ask for the time scale f* when the mutual information vanishes 

( 2 ) 


Equivalently, we ask for the time scale t* at which correlations between A and B vanish. 
When both subsystems A and B are the intervals {0 <)y < x < y + L and the perturbation 
is turned on at a; = 0 and time t_ = —we obtain the following analytical result for t* 3> /S 







2 sinh 


ttL 


(3) 


where (3 is the inverse temperature and = yl — 24/i^/c carries the information about the 
primary operator perturbation of conformal dimension h^. The parameter e represents a UV 
cut off for the local excitation, so that the excited state is localized around a region of size e 
of the operator insertion. This makes the energy of the perturbation hnite. In the 

limit h^/c -C 1, which is the relevant one to match the butterhy effect discussed in P [7], 
this reduces to 

C^/(i./?)4log(^). (4) 

where Sdensity = is the entropy density of the original thermal system. The log S behavior 
in Q is consistent with the fast scrambling conjecture PP. 

Given the bound Q, it should be possible to extract the same time scale from the 
condition of vanishing two sided 2-pt functions. We explicitly show this in appendix [C| 
conhrming the observation made in Shenker and Stanford P that both scales are controlled 
by the same physics. 

In the second part of this work, we derive the same time scale from bulk holographic 
considerations and hnd a perfect matching between both calculations. Our holographic model 

excitations are triggered by joining two semi-infinite lines and lead to local excitations in all sectors of a 
given CFT. Thus their behaviours differ from each other in integrable CFTs [131111117]. On the other hand, 
in large c CFTs, they behave similarly [miniiis] in that both results for 2d CFTs show logarithmic time 
evolution of entanglement entropy. 
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is based on the description of the local boundary perturbation in terms of some free falling 
particle satisfying an initial condition guaranteeing such particle carries the right amount of 
energy from the CFT stress tensor perspective. This is done by generalizing the model in 
PEI ED to the two sided BTZ black hole, based on the back reaction description of point 
particles as quotients of AdSs. Applying the holographic entanglement entropy [2S1 ESj to 
evaluate the entanglement entropy and mutual information in our set-up, leads to the same 
scrambling time (|^. At the same time, our setup and calculation may be useful for the 
interesting question regarding the dual CFT interpretation of a particle falling into a AdS 
black hole in future studies. 

Our local boundary perturbation includes a regularization parameter e describing its 
size. This parameter is holographically interpreted as the bulk position (distance from the 
boundary) from which the massive particle falls into the black hole. Our solution computes 
the back reacted geometry for any and approaches a localised shock-wave in the limit of 
large [23 E] • 

The paper is organized as follows: In section 2, we will analyze the time evolution of 
entanglement entropy in large c 2d CFTs at hnite temperature, which agrees perfectly with 
a previous gravity dual computation. In section 3, we study local perturbations in hnite 
temperature CFTs by employing the thermoheld double formalism. We compute the mutual 
information from entanglement entropies. In section 4, we compute the scrambling time for 
the mutual information. In section 5, we introduce our holographic model. In section 6, 
we present our holographic computations of mutual information in a two sided AdS^ black 
hole background with a local excitation. In appendix A, we explained the details of treat¬ 
ment of twist operators in the replica method computations of entanglement entropy in the 
thermoheld double formulation. In appendix B, we present some details of our holographic 
model. In appendix C we describe a computation of two point function in our model. 

Note Added: While hnishing our main computations, the work of Roberts and Stanford 
pn] appeared. The latter has a detailed account of two point functions in the presence of 
localised excitations over thermal states and briehy mentions the behaviour of the mutual 
information in the same set-up. Thus, it has some overlap with our results. In our paper, we 
literally evaluate the mutual information between the thermoheld double in both 2d large 
c CFTs and their gravity duals independently and show their results perfectly agree. Our 
gravity solutions explicitly have the regularization parameter e and our matching between 
gravity and CFT results holds while keeping this parameter small but non-zero. We would 
also like to mention that in the interesting recent paper [ST] by Maldacena, Shenker and 
Stanford, the fast scrambling behavior of the correlations functions has been interpreted in 
terms of chaos. 
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2 Single sided entropy 


To introduce our basic tools and fix the notation, we analyse the local perturbation to a 
thermal state in a single 2d CFT at finite temperature in this section. 

Consider a thermal state p/j locally perturbed by a primary operator ■^(0, —inserted at 
a: = 0 at time The time evolution of the resulting density matrix is given by 

p{t) = A/'e“*-^V(0, -tu,) pfi ^''*'(0, -t^) • (5) 


where H is the Hamiltonian of our system. 

Denote by pa = Tr^p(t) the reduced density matrix on a finite interval A with endpoint^ 
y + L satisfying y,L > 0. Its entanglement entropy Sa can be computed using the replica 
trick. We first compute the Renyi entropies 

( 6 ) 

1 — n 


The entanglement entropy is obtained by taking the limit Sa = lim„_,.i 

The trace of the reduced density matrix Trp^(t) requires the calculation of the normalised 

4-point function 


ry {'^{Xl,Xi)an{x2,X2)(Tn{x3,X3)'^^{x4,XA)) 


with the insertion points 


(7) 


xi =-ie, X2=y-t^-t_, X 3 = y + L - X 4 =+ie 

xi = +ie, X2 = y + t^ + t_, X3 = y + L + + t_, Xi = -ie. ( 8 ) 


The operator T stands for the product of the operators ipi in each of the i-th copies of the 
theorjj^ 

'h = • ^>2 • • • (9) 


and has conformal dimension = nhp, where is the conformal dimension of the original 
perturbation ip. Notice e is a parameter smearing the local operator perturbation and all the 
time evolution is carried by the twist operators cr^, which are initially inserted at both 
ends of the interval when cyclically gluing the different cylinder copies that give rise to the 
manifold C„. Finally, the conformal dimension = 2Hfj of the twist operators is 


H„ = 



( 10 ) 


We compute Q analogously to [2D] but with an additional composition of a map from 
the cylinder to the plane 

2tt 

w{x) = e '3 ( 11 ) 

^Notice that the perturbation is originally inserted outside of the interval A. 

^This correlator is formally computed in the cyclic or bifold CFT'^jZn- 
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to take care of the thermal nature of the original state, as well as the map 


z(w) = 


(wi -w)w 34 

Wl 3 (w - W 4 ) 


( 12 ) 


that brings the points wi —?■ 0, ^2 —t z, W 3 —)■ 1, tC 4 —)■ 00 . 

The transformation properties of primary operators determine the resulting trace to 


Trp^At) = 


. U f ^^23' 
-smh 






neuv V /3 
where we used the 2-pt function on the cylinder Ci 


{i){xx,Xi)i){x4,X4))ci = 

and introduced the canonical 4-point function 

1 4 / 1 . 


G{z, z) 


— sinh 

( 'XXi4\ 

TT 

K P / 


-4h,; 


G{z,z)= lim \Z 4 \*'^'^{'tp{z 4 ,Z 4 )an{z,z)an{lG)'^{^,^)) 

2:4—^■oo 


(13) 


(14) 


(15) 


dehned in terms of the cross-ratio 


2 ; = 


W 12 W 34 

W 13 W 24 ’ 


l-z = 


^ 141 x 23 

'U^13'?^24 


(16) 


We dehned Wij = Wi — Wj in all the above formulas and the same conventions hold for Xij 
and Zij. 

The Renyi entropies are computed by inserting (13) into (|^ 

1 


0 


^smh^)+- 
^TTEuv P n 


\og{\l-z\^^^G{z,z)), 


(17) 


where Suv is the UV cut off of the CFT i.e. the lattice spacing. The hrst term is the 
standard Renyi entropy for an interval L in a 2d CFT at hnite temperature T = 1/13] the 
second term captures the extra contribution due to the local operator insertion. In particular, 
the dependence on the conformal dimension of the local operator hy is encoded in G{z, z). 

In general, the extra contribution to the Renyi entropies requires the knowledge of the 
full four-point function G(z,z), i.e. the dynamical details of the particular 2d CFT under 
consideration. To make further progress, we consider the large c limitj^ Notice that in the 
limit n —)■ 1, the twist operators an, an become light 


Hfj/c ~ -/ ^ t 1 

^We already used the regularised twist operators so that euv is the standard UV cut-off. 
®We assume a class of CFTs allowing such a limit. 


( 18 ) 
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If h^/c remains fixed in the large c limit, the 4-pt function (15) becomes a 4-pt function 
involving two heavy and two light operators. This is precisely the set-up considered in 
mm to compute the dominant (vacuum) contribution to G{z, Using their results, one 
derives 1201 


_ c(n — 1) «V')(1 — 2 ;“’/>)(l — 

\ogG[z^z) ~---log 




+ 0{{n-lf) (19) 


where 

I 24^ , , 

= ( 20 ) 

encodes all the dependence on the conformal dimension of the local operator h^. Finally, 
we can compute the variation in the entanglement entropy due to the insertion of the local 
primary operator to be 





- z){l - z) 


( 21 ) 


where we subtracted the entanglement entropy of the interval L at hnite temperature T = 1//3 


‘S'thermal — g lo§ 


(——— sinh 
\7ieuv 



( 22 ) 


As explained in 121 , to extract a non-trivial contribution to the entanglement entropy in 
the CFT at hnite temperature we must take the smearing parameter e small but hnite. This 
still allows us to work with completely analytic formulas at order e. Then, in the small e//3 
limhj^ the cross-ratios are 


z ~ 1 -I- 


z ~ 1 


2TTie 


sinh ^ 


(3 sinh AlUhxUxM sinh 


2TTie 


sinh ^ 


13 sinh sinh 


+ 0{e^ 


+ 0(e' 


(23) 


Due to the non-trivial monodromy properties of G{z,z), we must carefully deal with the 
sign of the imaginary part of the cross-ratios [2D1 l3Uj^ Notice the imaginary part of z never 
changes sign, for -|-f^ > 0. Thus, we conclude ^ ~ 1 for all such times. On the other hand, 
the imaginary part of 2 ; does hip sign whenever t_ + E {y, y + L). Thus, we either have 

^We refer the readers to [321133] for the description of the regularisation of twist operators used to compute 
entanglement entropy. 

®We are assuming that y — t- — and y + L — t- — are larger than the smearing parameter e in units 
of/3. 

®We choose a reference phase to be consistent with causality and make entanglement entropies real and 
non-negative. 
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{z, z) —>• (1,1) for t- +tuj < y and t- +tuj > y + L or {z, z) —)■ (e^”, 1) for y < t + toj < y + L. 
Using these phases in (21), we reach our hrst important result 


ASa = 0 , < y and t- + > y + L 


A^A = ^log 
0 


13 sin Tia^j sinh 


sinh 


’ 

1 3 J 

\ 3 J 

726 

sinh 

(f) 



y<t_+t^<y + L. 


(24) 

Thus, there is no variation in the entanglement entropy Sa either till the perturbation reaches 
region A < y) or when it leaves region A > y+L). While the perturbation can 

causally be in region A, the variation in entanglement reaches a maximum at = y + ^ 


(ASa) 


max 



(3 sin na^ 
27ie 


tanh 


TT L 

2i3 


(25) 


In the high temperature limit (or large interval L), the increase in entanglement due to the 
perturbation equals 


ASa ^ X log 
0 


[3 sin 

27ie a^p 


[3 —^ 0 , — ^ 1 
e 


(26) 


In section |6.1| we will match this result with the gravity dual computation using the 
holographic entanglement entropy. Let us now proceed with the entanglement entropies that 
involve intervals in both CFTs. 


3 Two-sided entropies 


Entanglement entropy and mutual information in the thermoheld double state involving 
intervals in both CFTs were discussed in detail in [3Hj and their time evolution in [321 EZ] • 
Here we only briefly review this setup and extend by the insertion of a local operator to one 
of the CFTs. 

Consider two non-interacting 2d CFTs, CFTi and CFTr with isomorphic Hilbert spaces 
I-Ll and 'Hr. The thermoheld double (TFD) state is a particular entangled state in the total 
Hilbert space Htot = Hr® Hr. 

\'^3) = (27) 

where In) l,r G Hl,r are the eigenstates in each Hilbert space and Z(I3) is a partition function 
in total Hilbert space (and also in each Hilbert space). 


m =i: e 

n 

8 


— ^En 


(28) 


















If we have the total Hamiltonian Htot = Hi + this partition function is a partition 
function on a cylinder Ci with circumference (3 


n^m 


= Trtot 


(Hi^+Hr) 


(29) 


Tracing out the Hilbert space I-Lr from the pure state density matrix, we can get the thermal 
density matrix in CFTr with temperature (3 

Pi (30) 

A general time evolution of the TFD state is obtained by applying the evolution operator to 
both CFTs 

|T^(t_,4)) = (31) 

\JZ{I 3 ) n 

One can immediately check that setting t_ = = t, that can be seen as evolving with 

Hamiltonian Hr — Hr, leaves the TFD state invariant (the symmetry of the TFD state). 
On the other hand, setting t- = —t+ = t yields the time dependent state corresponding to 
evolution with Hr + Hr as in [32] • These two conhgurations should be kept in mind since 
we leave general and in our formulas so that our formalism can be used to extract the 
evolution after a local excitation with any of the two total Hamiltonian 

In the TFD formalism, we can also relate one-sided and two-sided correlators by the 
analytical continuation of t. For example, consider the following one-sided correlator 


{^fi\OL{xi,0)Ol{x2,t)\^l3) = ^^\n\LOL{xi,0)\m)L{m\LOl{x2,0)\n)L (32) 

n,m 

On the other hand, we can deform the two-sided correlator as follows 

0\OL{xi,O)OR{x2,t)\^ fi) = Xle"2(®"+^"*)+**(^"“®’-)(n|LC>L(xi,0)|m)L(n|HOij((C2,0)|m)R 

n,m 

= 2 {nlLOiixi, 0) |m)i,(m|ijC>Jj(a;2, 0) \n)R 

n,m 

(33) 


Therefore, the one-sided and two-sided correlators can be related through the analytical 
continuation t ^ t + 

{^p\OL{xi,0)OL{x2,m^) = {^p\OL{x,,0)OUx2,t-t^)\^p) (34) 

If the operators Or in CFTr are located at r = 0, the operators Or are located at r = |, 
or at the opposite side on the cylinder. We can express the correlators in the TFD state as 
the correlators on the cylinder Ci (see Fig. [^. 

fact our results also hold for evolution of the TFD state with Hr or Hr only and one can extract 
these formulas by setting t- = t and = 0 or = 0 and t+ = t respectively. 
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Figure 1: Our setup in the computation of the mutual information. We have two intervals 
A and B of size L 2 — Li = L m each CFT and a local operator inserted at time tw in the 
past. The operators are separated by distance 2e and in the CFT formulas we use Li = y 
and L 2 = y + L. 


Let us also mention a simple fact related to the symmetry of the TFD state. Namely, 
in the CFT, we compute the entanglement entropies as well as the mutual information in a 
state 

1^) = \^/J0) (35) 

with the TFD state I'lpy). Since — Hr leaves the TFD invariant, the above state is 
equivalent to 

1^) = . (36) 


Thus, the mutual information computed at = 0 in these two states has exactly the 

same functional dependence on This can be conhrmed from our explicit formulas for Ia-b 
in this section. 

Notice that in the previous single sided entanglement entropy calculations, we used trans¬ 
lation invariance to write the time dependence of the operator insertions in ([^ as a function 
of When computing two-sided observables, the same shift will be applied on the CFT 

time in the opposite boundary. This is consistent with the TFD path integral construction 
[3] and it also appears naturally in our holographic dual model as it can explicitly be seen 
in the embedding equations (160) and (161) appearing in appendix B. 


3.1 Semi-infinite intervals 

Before proceeding with hnite entangling regions, consider A and B to be semi-inhnite inter¬ 
vals X G [0, cx)). We want to clarify the difference between previous results for the second 
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{n = 2) Renyi mutual information in this setup [21] and our current mutual information 
[n = 1) discussion. In the large central charge c limit, and after the insertion of a local op- 
erator, the second Renyi entanglement entropy of the union S\(j^ grows linearly with time. 
Equivalently, the change in the second Renyi mutual information for semi-inhnite intervals 
decreases linearly with time [ 21 | 


A/( 2 ) ^ 


'nr 


(37) 


This holds for late times in the regime where 1 -C -C c. Below, we want to compare this 
behavior with a large c computation of the mutual information (n —)■ 1 ) with twist operators 
and for heavy local operators ~ C>(c) (as in [30]). 

To compute the entanglement entropy Saub between two semi-inhnite intervals A and B 
with starting point Li = ?/ > 0 on each boundary CFT, we must calculate 


Trp(^uB(^) 




(38) 


with the insertion points 


Xi = —ie, 
xi = +ie, 


X2 y to; t—, 


X2 = y + t^ + t-, 


^6 — y + i— — 

./9 

xe = y-i^ + t+ + t^, 


X 4 = +ie 


Xi = —ie. 


(39) 


Notice the edge of region B (the location xq) was shifted by 
the dependence on t+ is also through 


in accordance with (34) and 


We follow the same strategy as before : after mapping the cylinder to a plane hy w = e , 
it is the cross-ratio z, z on the plane that controls the relevant 4-pt function 


W 12 WM 

W1QW2A 

^12^64 

Wl%W2i 


~ 1 -|- 

~ 1 — 


2'Kie cosh ^ 

P sinh Ay-t--n) , pQgp Ay-i+^l 

2'Kie cosh At--t+) 

/3 sinh cosh 


+ 0{e^) , 
+ 0{e^). 


(40) 

(41) 


As before, the sign of the imaginary part of z hips for ti^ + t_ > y from positive to negative, 
whereas that of z is negative for all and t±. To the hrst order in e, we hnd {z, z) —)■ (1,1) 
for + < y and {z, z) —>■ 1) for toj + t- > y. Using the same method as for the single 

sided case, we obtain the entanglement entropy Saub for the semi-inhnite intervals to be 


Saub — 


flog 

flog 


cosh ^ 

TTeuv p \ 

(cosh ^ sinh cosh 

\7Teuv J ne 13 ^ /3 


(f_ +t^ <y) 

{y <t_+1^) 

(42) 
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where At = t_ — In particular, when t_ = t+ = 0 and is very large {y -C Sa\jb 
grows linearly with 


TTC C 

Smjb ~ “ S') + g 



^ /3 sin wap 

\ 7 reuv / 

47re ap 


( 43 ) 


We can hnd that the first term behaves like thermal entropy which is proportional to 1//9. 
Equivalently the mutual information decreases linearly with t^, but now with a coefficient 
proportional to the central charge c. As explained in m, this behavior is interpreted as the 
destruction of the entanglement between CFTi and CFTr and the broken "entanglement 
bond" reconnects between the subsystem AU B and its complement. The unit cost of this 
reconnection process is proportional to || (which is the entropy density Sdensity)- 


3.2 Mutual information for finite intervals 

In this section we compute the mutual information between hnite regions A and B in opposite 
boundaries in the TED state at large central charge. The setup is depicted on Fig. 

The mutual information is dehned as 

Ia-.b = Sa + Sb — Saub (44) 

where Savjb stands for the entanglement entropy of the union of the two intervals. Each of 
the three entropies is computed using the replica trick in terms of the correlators of the local 
operators T ([^ and twist fields inserted at the endpoints of the entangling regions. Since 
we already computed Sa, we focus on and Savjb- 


3 . 2.1 Sb 


The calculation of the entanglement entropy in the second GET/? is analogous to the one 
for Sa- It involves the same type of normalised correlation function 


AvpUt) 


(T(3:i, xi)an{x5, a^5)q~n(a^6, x^)) 


(45) 


but with the insertion points for the twist operators conveniently shifted by as reviewed 
above 


. ./3 ./9 

xi = -ze, X 5 = y + L + z- -t+ xq = y + z-- X 4 =+ze 

r -P 

xi =+ze, X5 = y + L - z-+t++ tuj, xq = y - z-+ X4 =-ze. 


(46) 


As stated earlier, is the shifted time being used in the right GET. To compute the 

4-pt correlator, we compose the map 0 with 


(wi - w){we - W4) 
{wi - we){w - W4) ' 


(47) 
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The corresponding cross-ratios equal 


27 rie sinh ^ 

P cosh cosh 


= ^5 ~ 1 


27 rie sinh ^ 

P cosh "(^+*++^-) cosh 


+ , 
+ 0 {e^). 


(48) 


Notice that the signs of the imaginary parts are the same for all and ty^. Thus, in the 
small e limit, {z, z) —?■ (1,1) for all t+. This reflects the intuition that the local perturbation 
turned on on the left CFT has no effect, at lowest order in e, in the quantum entanglement 
measured in the right CFT. Using the expansion in these cross-ratios, we derive 


glOg 


^-sinh^) VC 
TTSUV P / 


(49) 


Thus, quantum entanglement in the region B remains thermal for all f_|_ at lowest order in 
e, i.e. ASb = 0. 

3.2.2 Sa\jb 

The most interesting piece in the mutual information is Saub- Following [53], this requires 
the calculation of the 6-pt function 




{'Ipixi, Xi)an{x2, X2)an{x3, X3)an{x5, X5)an{xQ, ( 2 ) 4 , 2 : 4 )) 


{{'lp{Xi,Xi)'ljj^X4,X4:))ciT 
where the different insertion points correspond to the different interval endpoints 

2:1 = -ie, X2 = y - t- - tuj, X3 = y + L - t- - tyj, 2:4 = +ie 

xi = +ie, X2 = y + t_ X3 = y + L + t_ +tyj, X4 = -ie 

r 

2:5 = y + L + i--t+-ty,, Xq = y+ 

r -P - -P 

X 5 = y + L - i-+ t+ + ty„ XQ = y-i- + t+ + ty,. 


(50) 


(51) 


See Appendix for further comments on the ordering of the twist operators appearing in 

3 strategy as before, we comj 

w{x) = and z{w) = 


(50). Following the same strategy as before, we compose the two maps 

(wi - w)w3i 


Wi3{w - W4) ’ 


(52) 


and use the transformation properties of primary operators, to write the trace (50) as 
P . (7rL\ 


T^^Paub = 


|1 


U56 


an{z, Z)an{l, l)0-niz5, Z5)an(z6, ^^e)!^) 

(53) 
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where the cross-ratios (z,z) are given in (23), and Zi = z{wi). 

In the following, we discuss two different CFT channels: S and T-channel, where we 
compute this 6-pt function on the plane in the large c limit (see a detailed discussion in 
lED- The corrections to the particular channel choice are suppressed by e factors. We 
will explicitly see how these channels match the two different bulk geodesics determining 
the holographic entanglement entropy in our holographic discussions. The upshot is that S 
and T-channel correspond to the disconnected and connected geodesics for the holographic 
calculation of SavjBi respectively. 


S-channel : Let us introduce a resolution of the identity 

{'ip\(Tn{z, l)an{z5, Z5)an{ze, ZQ)\^jJ) = J2{'4’Wn{z, z)an{l, 1) |a) (a| an{z5, Z5)d-n(ze, ZQ)\^p) 

OL 

(54) 

where the sum runs over all possible intermediate states. 

Consider the hrst 4-pt function (-^|(Tn(^, 1) |q;). As we have seen in (23), the 

relevant limit corresponding to e —)■ 0, is either {z,z) —>■ (1,1) or {z,z) —)• (e^’^*,!). Thus, 
in either limit, the correlation function can be computed using the OPE of twist operators 


an{z,z)an{lA) ^^ + 0{{z-iy') r G (55) 

Ignoring the terms proportional to {z — 1) and focusing in the dominant contribution due to 
the identity operator, we reach the important conclusion that the summation over the entire 
set of intermediate steps is restricted to |a) = |'0) due to the orthogonality of 2-pt functions 
in any CFT. We stress that we could have reached the same conclusion in the limit of small 
L 2 — Li, but this is not required in our set-up. 

Thus, our 6-pt function can then be approximated by 

{'il!\an{z,z)an{l,l)'^n{z5,Z5)an{ze,ze)\'il:) ~ {'il!\an{z,z)an{l,l)\'4’) {'ll^\(^n{z5,h)an{zQ,Ze)\'ll!) 

(56) 

The hrst 4-pt function equals G{z,z) in (15), whereas the second 4-pt function factor will 


be proportional to the same function but evaluated at a different cross-ratio. To see this, 
consider the map taking zi —)■ 0 and Z 4 ^ oc 


z(x} = 


{zi - x){zg - Zj) 
(zi - Zg){x - Zi) 

This allows us to write the desired correlator as 

{^p\an{z5, Z5)an{ze, Ze)\^jJ) = |1 

Thus, the leading contribution in this channel is 


(57) 


(Tn{z5,Z^)(Tn{l,l)\'4’) ■ (58) 


TrPAuB - 


/S 


sinh ( ^ 

TlSuv V P / 


-SHa 


|1 - |1 - 2 : 51 ^^" G{z, z)G(z 5 , Z5) + ... 


(59) 


14 









where the dots stand for the contributions coming from the subleading terms in the OPE 
of the twist operators (55). Interestingly, since the cross-ratio equals z^, the cross-ratio 


determining S's, we reach the conclusion that 

Saub = Sa + Sb, 


and 


^A:B 


= 0 


(60) 


This channel reproduces the bulk expectation coming from geodesics joining points in the 
same boundary, leading to a vanishing mutual information. 


T-channel : We could also introduce the resolution of the identity as follows 

{^|J\an{z, z)an{l, l)(Tn{z5, Z5)an{Z(i, ZQ)\'ljj) = z)an{z6, Zq) Ict) (al Cr„(25, ^5)dn(l, 1)1-0) • 

a 

( 61 ) 

Notice the correlations involve twist operators inserted in different boundaries. Thus, we ex¬ 
pect this channel to reproduce the bulk contribution from geodesics connecting both bound¬ 
aries. Remember that in the small e limit, we already argued that 2:5 —?■ 1. Thus, we can use 
the same OPE argument as above to conclude that the dominant contribution comes from 
la) = {ip). By dehnition, this gives 

(01 crn(^5,^5)d„(l, 1 ) 1 - 0 ) =G{z5,Z5) (62) 

The remaining correlation is again proportional to the same function, but evaluated at a 
different cross-ratio. This is proved by considering the map 

X (zi - x)(z6 - Z4) 
z{x) = 7—^— --r . 

{zi - Zg){x - Z4) 

which allows us to derive 



\aniz,z)an{ZG,ZG)\^p) = 11 “ ^2 ^261 G{Z2, Z 2 ) 


\-4H^ 


(64) 


where Z 2 = z{w 2 ) = 2 ; as in (23) and Z 2 = z{z 2 )- Thus, after some manipulations we have 




— 

Tieuv V P , 


X 


1 — X 




1 — Z5\‘^^‘’\l — Z2 \'^^'’G{z2., Z2)G{z5, Z^) + .. 


I4n 


(65) 

where {x, x) are the cross-ratios computed out of the insertion points of the four twist oper¬ 
ators 

2 sinh^ ^ 


X = 




2 ^ 232^56 _ UI23WG6 _ 

2 ^ 252:36 W 25 W 3 G cosh -|- cosh 


= X , 


( 66 ) 


what allows us to write the dominant contribution from the T-channel as 

- 8 Ha 


TrpAuB - 


u f 


1 ~ ^ 2 ! '’G{z 2 tZ 2 )\'\- — Zg\ '"G{ZgtZg) + ..., (67) 
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where At = t_ — t+ and the cross-ratios 


^5 


= 1 _ _+ 0(,^, 

-^-' cosh ^(y+L-^^\ ^ 


/9 sinh 


2 me cosh 

= 1-1- P -u Oie^^ 

/9 sinh cosh ^iy+^+^++^^) 


/9 sinh ^ 


Z 2 = l + 


Z 2 = l- 


?!!!__ 

/? sinh cosh »<»-!-■■) ' 


( 68 ) 


27rze 


7r(t--t+) 


Now, nsing that at large 
have the identity [lUl IH] 


cosh- 5 — 

_^_u O(e^) 

/d sinh cosh ^ 

central charge and for two heavy and two light operators we 


\1-z\^’^Giz,z) ~ 


/z^2“ (1 — z°‘)z^2°' (1 — 2;“)\ 


as well as ( 68 ) we can extract the behavionr of Saub for any time regime, 
this carefnlly below assnming as before that 0 < y < y + L. 


(69) 

Let ns analyse 


It is clear that the monodromies of the correlator are determined depending on the relation 


of t- -|- tuj with y and y + L. From ( 68 ) the signs of the imaginary parts of z^ and Z 2 do 


not change with time and we have Z 5 ~ 1 and Z 2 — 1- On the other hand Z 2 — when 
t_ + > y and Z 5 ~ when > y + L. This gives ns three possible contribntions: 


•Saub ^ y log 

Q 2c 

Saub ^ y log 


/3 ^ fTTAG 

-cosh —— 

TTEuv \ P ^ 

(3 , ( vrAt' 


t_ +tuj <y 




cosh 


V p 


c , I j3 sin TTCK^ sinh 

log-— 

6 \ ire 


TT(t-+t^-y) , TT{t++t^-y) 

y y 


cosh ^ 


y<t_+t^<y + L 


(70) 


and for t_ -|- > 1 / -|- L we can rewrite onr trace as 


T^^Paub ^ 
(5 


13 (TTAt\ 

-cosh 


X 


TTSuV \ J 

[ttAG 




13 sin na^ sinh cosh ^ 


7r(t++tm-y) 


T^^uv 


cosh 


V /S 


TT ^ 

[3 sinTTcx^ 


cosh ^ 


sinh cosh ' 


- 2 Ha 


- 2 H„ 


TT e 


cosh^ 


(71) 
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The entanglement entropy Saub in this time regime can then be written as 


c (sinh cosh sinh cosh ^(^++yy -^) ' 

" 6 

2c. 


-log 


/3 / vrAt^ 

cosh 






cosh ^ 

t_ > y -\- L 


( 72 ) 


Notice that ( [7T| ) resembles the contributions from two different pieces. In the holographic 
part, these will be interpreted as the contributions from two bulk geodesics connecting points 
in opposite boundaries. 


3.3 The evolution of the Mutual information 

The evolution of the mutual information after turning on the local excitation can now be 
computed in the large central charge limit. At early times < ?/, the single sided 

entropies are thermal Sa — Sb = *S'thermai- Thus, the mutual information equals 


2c 


I sinh 


Ia-.B 3 >Og I jAi 


( 73 ) 


This is clearly hnite and depends on the Hamiltonian driving the evolution. If we use the 
bulk isometry then t_ = and At = 0, giving rise to a time independent 

mutual information, as it should. Notice that positivity of the mutual information in this 
case requires tiL/(3 > 1. Whereas for the Hb + Hr Hamiltonian, we recover the mutual 
information time decrease discussed in [SB]- 

These results can be understood using causality considerations : for < y, the 

perturbation did not enter into region A and could not possibly disturb the original thermal 


entanglement. Once the excitation reaches region A {y < t- + tuj < y + L), using (24), (49) 


and (70), the mutual information evolves as 


Ia-.b ^ 1% + X log 
0 


sinh Ay+L ^ —A) ^osh ^ 
cosh At++tB,-y) g- bL 


(74) 


Note how the dependence on the conformal dimension of the perturbation cancels, between 
the contributions in Sa and Saub, in this regime. It would be interesting to understand the 
mechanism behind this large c behaviour of the mutual information. 

In the last region t_+t^>y + L>y the mutual information equals 


Ia-.b — I%.B 



cosh At++t^-y) ginp .{t.+t^-y-L) ^^t+At^-y-L) ' 

cosh ^ cosh ^ 


( 75 ) 


17 































It is important to stress that when extracting the answer for the mutual information for 
various times as well as one has to maximise the mutual information between the S 
and the T channel answers so that it is always non-negative. 


4 Scrambling time 


Shenker and Stanford [0] dehned the scrambling time t* as the time scale at which the 
perturbation has destroyed all the preexistent correlations. In our notation, their condition 
reduces to setting = 0 and to study the vanishing of the mutual information 


-^A:s(C) - 0- 


Evaluating our previous results (73), (74) and (75) for = 0, we obtain 


Ia-.b ^ 


Ia-.b ^ 


Ia-.b ^ 


2c , ttL 
ylogsmh —, t^<y 


0 


(sinh ^)%inh^^ 






y <tuj < y + L 


(76) 

(77) 

(78) 


2c , L c ( (3 sin 7ra.0' 

— logsmh — - - log- 

3 p 6 \7re ^ 

'sinh sinh ' 


6 


t^> y + L 


(79) 


Notice the mutual information is a monotonically decreasing function of Thus, starting 


with a positive mutual information, i.e. vrL//S > 1, there is a single root t* where (76) holds. 
After that, by switching channels, the mutual information remains zero. 

The hrst question to answer is whether t* G (?/, ?/ -|- L) or whether > y + L. Clearly, 
the second condition can only hold if the mutual information is positive at the transition. 
This requirement gives rise to the constraint 


y I^ ^ lA:B{y + L + e) >0 


/3 / sin Trap' 


< 


TTC 




sinh^ ^ 
cosh ^ 


(80) 


where we already used e -C /3, as in our previous CFT analysi^^ Because of working in 
this region of parameter space (e -C /3), we conclude that only small perturbations (ap —)■ 1) 
allow scrambling time scales t* > y + L. Since the function sin(7rap)/ap is monotonically 
decreasing in ap (or increasing in hp), the smaller the perturbation is, the easier it is to fulhl 


condition (80) for generic values of L/jl. Since this is the regime considered in [6], we will 
study the scrambling time under these circumstances^ 


e /? breaks down, then condition (801 is modified. 

^^One can equally study the conditions under which t* G (y, y + L). These generically require heavier 


perturbations. We do not understand this regime, which appears precisely when the mutual information (741 
is Up independent. 
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For > y + L and At = 0, the mutual information (75) becomes 


^ 77 log 


p \'Ke j 


-2 


6 cosh y) y l) 


Trjtuj—y-L) 

a 


( 81 ) 


This vanishes when 

/ 2'Ka^esmh^ 

I /Ssin(7ra^) 


= sinh^ ~ cosh ^"^ 


/S 


(3 


1 - 


tanh ^ 


tanh 


/3 


(82) 


Notice that — y > L guarantees the positivity of the left hand side. Condition (82) gives 
rise to a quadratic equation in sinh^ ^Trp^-y) ^ shown that there is a unique consistent 

root, in agreement with our previous arguments. In the limit t*//5 3> 1, this root reduces 
the scrambling time f* to 


C-!'+ 2 - ^ 


^ /ysh^\ ^ 

^TT \7re J TT y p y 


(83) 


Due to the non-compactness of the 2d CFT, no recurrences were expected to be seen in our 
calculation. Working in the small h^/c limit, as required by our analysis, then 


/3 sirnia.^ 

7 V€ (y^p 


E.. 




5, 


where S, 


density 


density 


TTC 

ibd’ 


(84) 


and is the total energy of our local excitation given by integrating the energy 

density as in ng. In this limit, the scrambling time reduces to 


1_ f ^density P i _ _ (r\ _ ■ i ^ 


2 4R, 


-) 


-I- — log 2 sinh — . 


TT 






(85) 


The log S dependence is indeed consistent with the original scrambling conjecture ISIS]. 


5 Holographic description 

In this section we compute the mutual information using the AdSa gravity dual of the 2d 
CFT set-up discussed in the previous section. The starting thermal state pp in the quantum 
theory is described in the semiclassical approximation of the correspondence by a black hole 
in the gravity side, the BTZ black hole [H]. Observables involving a single Hilbert space are 
described by the metric outside of the event horizon, whereas observables depending on both 
Hilbert spaces require the Kruskal-like extension of the BTZ black hole, as described in |3]. 

To compute the time evolution of the holographic entanglement entropy, we approximate 
the local CFT perturbation with conformal dimension A(= 2h^) by a bulk free falling massive 
point particle with mass m = A/R [TS]. Next, we compute its back-reaction on the BTZ 
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background in Kruskal coordinates using the coordinate transformation in Finally, 

we compute the entanglement entropy and the mutual information using the holographic 
prescription [ 2 H 1121 ] • 

In three dimensions, the back-reacted metric of a point particle at r = 0 in global coor¬ 
dinates is known 

= -(r‘^ + - /i) o ^ -^ r'^dip'^ , (86) 

where the mass of the point particle is related to p by /i = SGNR^m = ^ jg 

radius of AdSa. Depending on the mass of the particle, the background describes a conical 
singularity or a BTZ black hole. 

Now, the holographic entanglement entropy is given by = | times the length 
of geodesic 7 which connects the boundaries of the subsystem A for which we dehne the 
entanglement entropy Sa- In the above metric, the entanglement entropy of the boundary 
region A with endpoints and is |IH] 

Sa = ^ log 
6 

where a = carries the information on the perturbation, as in the CFT 

discussion, Atoo = and A(poo = satishes 0 < |A(poo| < tt. 

Mapping the static r = 0 geodesic to one starting at some distance e from the boundary 

and falling into the horizon afterwards can approximate the local perturbation turned on 

in the boundary theory. This is precisely the approach followed in [25 to describe the time 
dependent evolution of entanglement entropy in the bulk for locally perturbed thermal states. 
To describe the evolution across the horizon, which is required to study two sided correlation 
functions, one must use Kruskal coordinates. This is one of the tasks we will undertake in 
this section. 


2 rW . j-m cos (|Aroo|a) - cos (|Av?oo|a) 

a? 


(87) 


5.1 Free falling particle in Kruskal coordinates 

The geodesic of a free falling particle in the AdS-Schwarzschild patch of the BTZ black hole 




riA 

- (1 - Mz^) dtl + _ 2 + 


Mz^ 


was already computed in [ 21 ] to be 

t_=f, 9 = 0, 1 - Mz^ = (1 - Me^) cosh “2 


6 ^ 6 2ti 


M(T + tA) ■ 


( 88 ) 


(89) 


The only addition in the expression above is the shift r —)■ r -|- to account for the initial 
boundary condition z{—ti^) = e guaranteeing the particle’s energy 

m R 


E = 


Vl - Me2, 


(90) 
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matches the energy of the CFT perturbation in the small e limilj^ We extend this result to 
the entire eternal black hole by working in Kruskal coordinates. 


One way to achieve this goal is to map the global AdSa description (86) to Kruskal 


coordinates. A second one is to solve the geodesic equation directly in Kruskal coordinates. 
We check below that, as expected, both approaches agree. 

Free falling particle in Kruskal coordinates : The Kruskal extension of the BTZ 

, -4dT2 + AdX^ + (1 - T2 + d(j)^ 


metric (88) is given by 


2 2 + (—1 + 2 


ds^ = W 


(1 + UvY 


= R^ 


{1 + T^-x‘^y 


(91) 


where u = T — X gM, n = T + X gM with their range satisfying — 1 < mu < 1 and 
0 ~ 0 + 47r^//5. The conformal boundary, horizons and singularities are at uv = —1, nn = 0 
and uv = 1, respectively, with the left and right Kruskal regions dehned by 


Left: R- = {0 < u, —1 < nn < 0} 

Right: R+ = {m < 0, —1 < mu < 0} 

Both coordinate systems are related to the AdS-Schwarzschild patches via 


(92) 


V zh + z 


T = ±, 


1 — y/Mz . , / r— 

\ -^^smhfvMtqr 

^ 1 + VMz ^ ^ 


n = =Fi 




\ ZH-z ^_,^l,^ 
Zh + Z 


1 — y/Mz , / r— 

\ -^^coshfvMti 

\ 1 + y/Mz ^ ^ 


(93) 




coshl vM 


0 ) 


Using these maps ( p^ , we can rewrite the geodesic (89) in the parametric form 
A'(r) = 

T(f) = 


1 + ^/l — (1 — Me^) cosh ^ (^^/M (r + U 


(94) 


coshl V M 


d) 


1 + — (1 — Me^) cosh ^ (^y/M{T + 


The initial condition (t_,z) = is mapped to 


(TthXo) = ^ (-sinh (Vmc) ,-cosh 


(95) 


the absence of matter fields and working in the large c limit, as we do, the satisfaction of this condition 
guarantees that our holographic model should capture the bulk description of the identity conformal block, 
which in 2d includes the stress tensor. 
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This allows us to determine = t^(To,Xo). Similarly, r = t(T,X) can be determined from 
(93). Altogether, we can solve for T{X) as 


T(X) = 


sinh 


± 


X + 


cosh (a/m\ Me^ 


Vl - Me2 


1 -Me2 


Vl - Me2 \ 

Proceeding in a similar way, we can obtain the geodesic v = v{u) that reduces to 

ttiU — 1 


with 


Ui = 


1 - 


v{u) = 


U -\- 0,2 

1 - UqVq 

CL2 — 




(96) 


(97) 


(98) 


2'Uo VI - M€2’ 2'tio VI - A7V ' 

Checking equations of motion : Consider the relativistic action for a particle of mass 
m moving in the background metric (91) at constant 0. Working in the gauge where the 

(99) 


parameter along the curve equals u, this effective action reduces to 

^/v'du 


Its equation of motion 
has general solution 


S = -2mR , 

J 1 + uv 

v"{uv + 1) — 2v'{uv' — u) = 0 , 


( 100 ) 


, , C 2 + (C.+C|)'t. „ , Cl 

= - --T T c i r -- - ” <“' = (TTc^s 


( 101 ) 


This agrees with the geodesic solution (|97|) if the integration constants are matched as 

2v[ 


Vo{uQVo + lf ,.2 2 x/Mt. 


^2 = 


= -Vl - 


1 - UqVq 

In T, X coordinates these extended geodesics have two branches which meet at the point 


( 102 ) 


^m. — 


y/Me — cosh (^'/M 


Vl - Me2 

They cross the future and past horizons at 


T = 

-*■ m. 


sinh (y^t^ 
\/l-Me2' 


Xh± = 

and hit the past and future singularities at 

Ao± = ±sinh fv^MT 


(103) 


(104) 


(105) 

All these features can be seen in the sample geodesic plotted on the T — X Kruskal diagram 

(Fig§. 
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Figure 2: Plot shows our time-like geodesic on Kruskal diagram. The red part is given by 


(94) and the full geodesic (blue) by (96). Plot for M = 10, e = 0.01 and = 0.25 


5.2 Back-reacted metric 

The back reaction of the free falling particle in Kruskal coordinates can be obtained by 


following P! and rewriting the metric ( [M| in Kruskal coordinates, but taking into the 
account the specihc initial conditions discussed above. To solve this problem when = 0, 
in [21], we considered a boost in the plane Xi — X 3 with rapidity A 2 = A 2 (M, e) . In order 
to introduce the further parameter we will consider a preliminary boost in Xq — X 3 , 
since this corresponds to the natural boost action on the light-like coordinates u — v which 
captures the blue-shift of energy near the horizon stressed in [H|. In practice we then apply 
two particular boosts into the embedding coordinates of AdS^ such that the identihcation 
between global and Kruskal coordinates becomes 


sin T = cosh Ai Xq -|- sinh Ai X 3 
+ e~^^v 


= R- 


1 + uv 


cos r = cosh A 2 Xi — sinh A 2 (sinh Ai Xq -|- cosh Ai X 3 ) 
i?cosh A 2(1 — uv) 


cosh 0 — tanh A 2 


e 


uv 


1 + uv 

r sin ip = X 2 = R -sinh (j) 

1 + uv 

r cos p = — sinh A 2 Xi -|- cosh A 2 (sinh Ai Xq -|- cosh Ai X 3 ) 
R cosh A 2(1 — uv) (e^^u — e~^^v 


1 + uv 


uv 


— tanh A 2 cosh 0 


(106) 
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Solving for r = r{u,v,(j)) 


r = 


R{1 — uv) cosh A 2 


sinh^ 0 — e~^^v , , , 

A - 5 -h-tanhA2Cosh0 (107) 

\ cosh A 2 \ 1 — / 


1 + uv 

we can determine Ai and A 2 reqniring that the location of the static particle in global AdSs, 
r = 0, gets mapped into the free falling geodesic v{u) in (97). This hxes both boost param¬ 
eters to be 

Xi = , tanh A 2 = Vl — Me^ . (108) 

These boost parameters determine the explicit map between global AdSa and a free falling 
particle in Kruskal coordinates : 


r = 


R 

1 — UV 


VMe 

1 + uv 

\ 


and also 


Me^ sinh^ 0 -|- 


tan r = a/M e- 


g'v /DAttjj _ g—\/ A4tij 

1 — UV 




Vl — Me^ cosh (j) 


(109) 


1 — UV 


tan (p = a/m e- 


cosh (f) — \/l — Me^ 
sinh(^ 


\ — UV 


)%/ iVlt(jj g — y/Ki 

1 — UV 


— a/1 — Me^ cosh cf) 


( 110 ) 

( 111 ) 


Using this map, we can compute the exact back-reacted metric corresponding to a free falling 
particle in the eternal black hole satisfying the initial condition (uq, Vq). 

Our analysis is valid for any value of toj. This allows us to compare with some approaches 
in the literature where the back reaction of the local perturbation in the CFT was approx¬ 
imated by a shock-wave, i.e. a BTZ spacetime in the presence of some non-trivial stress 
tensor localised at the horizon. Since our approach in 3d was based on computing the ex¬ 
plicit backreaction of some point particle moving in some geodesic into the BTZ geometry, 
we can study the limit S> 1 in our geodesic analysis in subsection 5.1. In particular, £g- 
ure illustrates how our particle geodesic approaches a null geodesic on the horizon for such 
large Thus, our back-reacted metric in this particular limit should indeed correspond to a 
shock-wave propagating in the BTZ background as originally described in |H]. The advantage 
of the shock-wave description is that it also applies in higher dimensions, whereas our hnite 
results show the agreement between CFT and bulk computations also hold beyond this 
regime. 


6 Bulk mutual information 

The mutual information Ia-.b between regions A and B in the left and right boundaries, 
respectively 

Ia:B = Sa + Sb — Saub , ( 112 ) 
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can now be compnted by applying (87) to the three different bnlk geodesics providing the 
relevant minimal snrface compnting entanglement entropy in the bnlk. All we need to know 
are the locations of the endpoints in the limits of small e and Zqo that we will insert into (87). 
To keep the gravity formnlas compact, the endpoints of the intervals will be denoted by Lj, 
i G {1,2}, where 0 < Li < L 2 . To compare with the CFT formnlas one can snbstitnte 
Li = y and L 2 = y + L. 


6.1 Geodesic in the left boundary 

The two endpoints of the entanglement region A in the left bonndary are {t_,zi,6i) = 
{t_, Zoo, Li) and {t_,Z2,9i) = {t_, z^o, L 2 ). It is convenient to compnte their image in global 
AdSa nsing the asymptotic maps to the right and left regions (see Appendix]^. Proceeding 
this way, their radial coordinates in global AdSs satisfy 


^(1)^(2) ^ ^^ D 1 D 2 (113) 

\MeZooJ 

where 

Di = \ cosh VmL i — cosh VM{t-+ ti^)\ i = 1,2 (114) 

whereas the other coordinates are 


tan 

tan 


VMe 

VMc 


sinh 1 



1 

cosh 1 


— cosh 1 

(^\/ M (t_ + 


sinh 1 

[VMLi] 

1 

cosh 1 

M (t_ + 

1 — cosh 1 

{VMLi) 


(115) 

(116) 


with i = 1,2. 

The length of the geodesics depends on the valne of the time argnment t- + tw When 
+ tw < Li < L 2 , then the bonndary points eqnal 


T 


~ Tm, 


sinh 


M (t + toj)) 


D, 



These determine the coordinate intervals to be 


( 117 ) 


lArl ~ 


|A(p| ~ 


^Me 

D1D2 

^Me 

D1D2 


D 2 — Di \ sinh + t^), 

Di sinh L 2 — D 2 sinh ^/MLi 


(118) 


Dne to the identity 


DiDsdAipl^ - |Ard) = 4Me^sinh' 


71AL 


(119) 
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the geodesic length is 


~ log 


2rd)r(2) 

cos (a| Ar|) — cos (a| A(y9|) 
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~ log 


(l)r.(2) 






-(|Ay,p-|Arp 


~ 2 log 


/S 


sinh — 

.TTZoo P , 


( 120 ) 


This reproduces the thermal entanglement entropy (22) computed in the CFT in the same 
time interval once both UV cut-offs are identified euv = ^oo- 
Similarly, when > L 2 > Li, the boundary points equal 


j{^) ~ TT — a/m e- 


sinh -|- tij)) 




~ a/M( 


sinh (^y/MLi^ 


D, 


( 121 ) 


These are different from (117), but give rise to the same intervals (118). Thus, the length of 
the bulk geodesic joining them equals (120). This matches our CFT again. 

sinh \^MLi 


Finally, when L 2 > t_ -f > -^ 1 , the boundary points are 




^(2) ^ y^^siiffiA/M(f_^t^ ^ 
D 2 




(2) ^ ^ _ a/Mc 


Di 

sinh '/ML 2 
~D2 


( 122 ) 


From them we can easily get the absolute values of the intervals 

^/Me 


|Ar| ~ TT 
|A(p| ~ TT 


DID 2 

y/Me 

D 1 D 2 


{Di + D 2 ) sinh a/m( t_ +t^), 
(^Di sinh a/mL 2 -|- D 2 sinh a/m L 


(123) 


Notice that in the small e limit we are working on, |Ar| and |A(p| are close to each other 
^/Me r 


5 = |A(y9| — I Ar| = 


(Hi -|- D 2 ) sinh a/m {t- + t^) — Di sinh v/VfL 2 — D 2 sinh v^MLi 

(124) 


D 1 D 2 

This allows us to write the length of the bulk geodesic between these two boundary points 


as 


~ log 


~ log 
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~ log 
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-sinh ——- - -^^ 

TTZoo p ) TTe a sinh 


(125) 


where AL = L 2 — Li. This also perfectly matches our CFT result (24) after employing the 
Ryu-Takayanagi formula. 


26 







































6.2 Geodesic in the right boundary 

The two endpoints of the entanglement region B in the right bonndary are (t+,Zi,0i) = 
(t+, Zoo,-f^i) and (t+,Z 2 , 6 'i) = (t+, Zqo,-^" 2 )- Their radial coordinates satisfy 


^( 1 )^( 2 ) ^ 


f—) 

V MeZoo / 


D 1 D 2 


where 


Dj = I cosh+ cosh \/M(f++ i = l ,2 
whereas the other coordinates are 

sinh (^\/M (t+ + 


tanr*-*^ ~ —\/Me 


, tan(y9^*^ ~ 


,_ smhi^/ML^ 

-VMe -^ 


A " D, 

In this case, no matter what the value of is, the boundary points are identihed as 

sinh 


(126) 

(127) 
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These give rise to the intervals 

y/Me 
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D 1 D 2 
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Di- D 2 \ sinh a/M( t+ + t ^), 
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(130) 


Using the identity 


- I Arp) = 4MeAinh' 
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the geodesic length equals 

2 rd)rU) cos (a| Ar|) — cos (a| A(p|) 


~ log 
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(132) 


This reproduces the well-known thermal answer obtained in the CFT 


A ^ g log 
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vrzo 


7rAL\ 


= S, 


thermal 5 


(133) 


which also agrees with the CFT expression for Sb in (49). 


27 























6.3 Geodesics across the horizon and Mutual Information 


We want to compute the geodesic length between two opposite boundary points located 
at the same space like location but with different time labels We will describe the 
calculation once and apply it to the two cases of interest afterwards. The product of the 
radial coordinates equals 

R 


^(1)^(2) ^ 


f—y 

V Mezoo J 


D1D2 


(134) 


where 


Di = I cosh — cosh vyM(t_ + t^)| (135) 

D 2 = I cosh v^MLj + cosh vyM(t++ t^)| (136) 

where L* labels the space like location in both boundaries, i.e. either Li or L 2 . The other 
coordinates for the left boundary point are 

sinh v^M(t_ + 


tanr*-^^ ~ v^Me 
tan ~ y/Me 


cosh v^M Li — cosh ^/M (t_ + t^j) 
sinh y/MLi 


cosh v^M (t_ + toj) — cosh v^M Li 
At early times, L* > these are given by 

sinh v^M(t_ + tui) 
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whereas at late times. 
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VMe 


sinh v^M (t_ + 


cosh {t- + to;) — cosh Li 
sinh v^MLi 


cosh v^M [t- + t^) — cosh v^M Li 
The remaining coordinates for the right boundary point are 

sinh (t+ + tai) 
cosh v^M Li + cosh v^M(t+ + 
sinh Li 

cosh v^M(t+ + tu,)+ cosh v^ML* 

In this case, they are always given by 
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(141) 
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(144) 


7-(2) ^ 


sinh y/M {t+ + tu 




( 2 ) ^ ^ 
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(145) 

Do 
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(146) 



Let us compute the length of the geodesic in the early time regime Li > In this case, 
the interval differences are 


\At\ = — 

|A(p| = — 


VMe 

D 1 D 2 

VMe 

D 1 D 2 


D 2 sinh Vm (f_ + tuj) + Di sinh (f+ + toj) 
D 2 sinh y/MLi — Di sinh y/MLi 


(147) 


Plugging this into the geodesic length (|87|), we obtain 
~ log 
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In the late time regime, the interval differences equal 

y/Me 


(148) 


|Ar| ~ TT ■ 
|A(p| ~ TT 


DID 2 

y/Me 


D 1 D 2 

Since they are very close, we have 
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(^Di sinh y/M Li + D 2 sinh y/M Li 


(149) 
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(152) 
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This allows to write the geodesic length as 

cos (a|Ar|) — cos (a|A(p|) 
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(153) 


where At = t- —t+. These geodesics can now be used to compute the entanglement entropy 
of the union Savjb- 

In particular we will be interested in large > L 2 > Li when Sa = Sb = •S'thermai- In 
this case, there is a competition between the two geodesics connecting points in opposite 
boundaries and the geodesics connecting points in the same boundary giving rise to 2Sb- 
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The length of the new geodesics is 


Li ~ log 


~ log 


- ■ ,sinhcosh 


(3 cosh \ ^ sin na ' 


71Z„ 


Tie a 


cosh^ 


' /3 cosh ^ y /3 sin 7 ra sinh cosh ^ 
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7re a 


cosh 


( 154 ) 

( 155 ) 


where again At = t_ — t+. 

Summarizing, the holographic entanglement entropy of the union of two intervals on the left 
and the right boundary is given by 


SauB - F (-^7 + -^ 7 ) ’ 


(156) 


which matches with the CFT result (72). 


Finally from the above holographic results, we obtain the holographic mutual information 
= Sa + Sb — Saub and this again reproduces the CFT result (75) perfectly. As a 


B 


IA: 

consequence the scrambling time derived in the CFT (|83|) also holds as a result in gravity 


replacing —>• a, as stressed below (87). 
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A Twist operators in the TFD state 

In this appendix we explain the proper ordering prescription of twist operator insertions 
when computing the two sided entanglement entropy Saub (or mutual information Ia-.b) in 
the thermofield double state. 

When computing Saub, with regions A and B belonging to opposite boundaries, the 
replica trick instructs us to insert twist operators (T„, on each boundary. The order of 
these insertions is important, because different orderings can give rise to different replica 
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geometries. There are two kinds of insertion orders. One, where both bonndaries have the 
same order, snch as 

(^ 71 ) ^n)Li ^n}R • 

In this case, going aronnd the replica n-sheeted cylinder by passing throngh the cnts between 
the twist operators on each bonndary, one retnrns to the starting point after going aronnd 
n cylinders (see the right pictnre in Figj^. There exists a second insertion order in which 
both bonndaries have opposite twist operator orders, snch as 

(^ 71 ) ^n)Li {.^ni ^n)R ■ 

The same operation as above returns to the same point after going once around a standard 
cylinder (see the left picture in Figj^. 
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Figure 3: In the left picture (the second choice (158)), one returns to the starting point 


after going around two standard cylinders. In the right picture (the hrst choice (157)), this 
operation involves going around n standard cylinders. 


The question is what the right order prescription is when we consider the TFD state. To 
answer this question, we provide two CFTs and one holographic bulk arguments. 

It is convenient to remember the TFD state can be constructed by Euclidean time evo¬ 
lution from the pure state Yl,n \^)L\n)R with Euclidean time ±/S/4 for each boundary. To 
compute the density matrix paus, one glues partially by imposing boundary conditions to 
each sheet (half cylinder) . Finally, to construct Trpdu^, one glues the cylinders through 
their cuts on the cylinders, giving rise to a partition function on the n-sheeted cylinder. If 
one goes around this n-sheeted cylinder through the cuts, one returns to the starting point 
after going once around a standard cylinder because one should obtain the original state in 


the limit /3 —?• 0. Thus, the second insertion order (158) for twist operators is selected in this 
way. 
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A further argument to confirm this choice is as follows. The n-sheeted replica geometry 
is constructed by gluing single-sheets together along the cuts. If one enters the cuts from the 
negative region, one emerges in the upper sheet, whereas if one enters from positive region, 
one ends in the lower sheet. From this gluing condition, one obtains the n-sheeted replica 
geometry corresponding to the second type of insertion order of twist operators (see Fig. |^. 



Figure 4: In the replica geometry which we obtain by gluing n cylinders along the cuts, we 
can go around 2 standard cylinders, for example, Ath and {i -|- l)-th cylinders (dotted line). 

Our hnal argument is holographic. In the bulk, one is instructed to consider geodesics 
connecting the different edges of the subsystems. Twist operators are inserted in these edges. 
For geodesics connecting edges on the same boundary, the inserted twist operators must 
have the same monodromy (or charge) properties, whereas for geodesics connecting edges 
in opposite boundaries, twist operators must have opposite monodromy (or charge). For 
example, if cr„ is on an edge, an must appear at the other edge. This consistency condition 
chooses the same twist operator insertion order as in the CFT side. 

B Details of the holographic model 

Here we collect several useful formulas and conventions that we used in the part with holo¬ 
graphic computations. 

The relation between Kruskal and AdS-Schwarzschild coordinates can be obtained by 
referring both descriptions to the where AdSs becomes the quadratic surface 

- - Xl + Xl + Xl = -R^ (159) 
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in which we have 
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2X 
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y/Mz 

This formulas hx our conventions for the appropriate signs on the gravity side. 


(160) 


In order to derive lengths of the geodesic in the back-reacted metric we only need to 
know the asymptotic form of the map from the left and the right wedge to AdS's in global 
coordinates. This is just the map that takes the trajectory of our massive point particle in 
Kruskal coordinates to the r = 0 particle in global AdS^. In the left and right exteriors the 
map becomes 


+ r’2 sin T = 
V+ r'^ cos r = 
r sin ip = 
r cos (f = 


rJi - Mzl^ f — 

' ■ nX -■ + tR) 

y/Mzoc ^ 

(^cosh [^/Me) T Vl- Mz^ cosh , 

±0^1 — Mzl^ cosh + t^J)) — a/I — Me2 cosh [yjAB^ . 


R 

Mezoo 

R 


R 


Mez, 


(161) 


The points in the main text are extracted to the hrst order in e. 

Note that the map from the right wedge also depends on e as well as 


C Two sided 2-pt functions 


Given the bound ([^ that the mutual information provides on the amount of correlation, it 
is natural to study the two sided two-point correlation function in the TFD after a local 
perturbation is turned on in one of the boundaries. This is the calculation described in this 
appendix. 

Given a local primary probe operator Oh{x,x) of conformal dimension h and a pertur¬ 
bation described by a different primary the appropriate normalised 4-point two sided 
correlation is 


^ _ {Oh^ {Xl, Xi)Oh{x2, X 2 )Ohix 3 , X3)Oh^ {X4, 3 : 4 )) 
^ {Oh^{xi,Xi)Oh^{x4,X4)} 


(162) 
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where the insertion points are 


Xi = —le 


Xi = le 


X2 — Li — t_ — X 3 — L 2 


./9 
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X2 — Li + X3 — L2 + — i — 


Xi = —%€ 


(163) 


Assnming h/c -C 1 and keeping h^/c hxed, the nnmerator is again a 4-point fnnction involv¬ 
ing two heavy and two light operators. Using the large central charge resnlts of ra we can 
write the correlator as 

-Ah 
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13 . , vra;23 

-sinh —— 


-KZr, 


1 — G{z, z) 
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In the limit of small e/(3, the cross-ratios rednce to 
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(164) 


(165) 
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(167) 


where AL = L 2 — Li and At = t_ —1+. As in onr main text discnssions, the key observation 
is that the imaginary part of z depends on the sign of Li — t_ — For Li > t_ +we 
have (z, z) ( 1 , 1 ) and 

-Ah 


a- 




\2 


1 ( 2tiAL 271 At\ 

- 1 cosh -;;-h cosh 


/3 
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whereas for Li < t_ + {z, z) ~ (e^”, 1 ) and 
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1 ( , 27rAL , 27rAt\ 
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sinh A ^ 


Tie a 


cosh UAl+ap 


(169) 

It can be checked this resnlt precisely matches the gravity compntation where the two-point 
fnnction is given by the length of a geodesic (153) between two-bonndaries in onr back-reacted 
metric. 

One can read off the scrambling time scale from the two-point correlators. For example, 

3 r large 

(3 sin(7rQ;) 


setting Li = L 2 = 0 and t_ =t+ = 0, the correlation (169) for large is given by 

-Ah 


04- 
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27lZn 


exp 


dvrh / (3 
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^Tie a ) ^ 

This reprodnces the dependence on the pertnrbation for the scrambling time derived in the 
mntnal information analysis. 
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